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ABSTRACT 

Aims. Exoplanet-host stars (EHS) are known to present superficial chemical abundances different from those of stars without any 
detected planet (NEHS). EHS are, on the average, overmetallic compared to the Sun. The observations also show that, for cool stars, 
lithium is more depleted in EHS than in NEHS. The aim of this paper is to obtain constraints on possible models able to explain this 
difference, in the framework of overmetallic models compared to models with solar abundances. 

Methods. We have computed main sequence stellar models with various masses and metallicities. The results show different behaviour 
for the lithium destruction according to those parameters. We compare these results to the spectroscopic observations of lithium. 
Results. Our models show that the observed lithium differences between EHS and NEHS are not directly due to the overmetallicity 
of the EHS: some extra mixing is needed below the convective zones. We discuss possible explanations for the needed extra mixing, 
in particular an increase of the mixing efficiency associated with the development of shear instabilities below the convective zone, 
triggered by angular momentum transfer due to the planetary migration. 



1. Introduction 

More than one year after the launching of the satellite CoRoT 
and at the time we receive the first data, the understanding of 
the formation of the planetary systems becomes crucial. In the 
present state-of-the-art, the differences in the protostellar clouds 
and in the internal structures between exoplanet-host stars (EHS) 
and stars without an y detected planet (NEHS) are not well- 
know n. Many studies (iGonzalezll 1998t ISantos et aLll2000i |200TI 
|2003|) have shown that EHS are on average more metallic than 
NEHS. Two scenarii concerning the formation of planetary sys- 
tems have been proposed to account for this behavior: an original 
overmetallic protostellar cloud or a process of accretion of over- 
metallic matter (planetesimals, asteroids, comets, etc.) onto the 
star. This second scenario is now ruled out for several reasons 
listed below. 

Prec ise spectroscopic ground-based observations o f 
lithium (Israelian et al.ll200l lRvanl l2000t llsraelian et al Il200 
IChen & Z hao 200^ and beryllium jSantos et al.ll2002l 12004" 
in EHS and NEHS provide new constraints for stellar modelling. 
The study of the light element abundances may give us addi- 
tional information about the rotation and angular momentum 
history of stars, particularly for the EHS. A different angular 
momentum history, and consequently different light-element 
abundances can be expected for EHS and NEHS, due to the 
presence of a more massive protostellar disc and possible 
accre tion of planetary-mass bodies onto the star dSiess & Livid 
[19991) . 



Israelian et al.l(l2004l) show that a significant difference in the 



lithium abundance between EHS and NEHS in the range of ef- 
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fective temperature 5600-5850 K does exist. In this range of ef- 
fective temperatures, EHS are more lithium underabundant than 
NEHS. 

In this paper, we compute various stellar models to com- 
pare the theoretical lithium abundances to those determined from 
the spectroscopic observations. In the first part, we present the 
metallic peculiarities of EHS. In the second part, the evolution 
of various types of models calculated with the Toulouse-Geneva 
Evolution CodeJTGEC) are presented and c ompared to the ob- 
served stars of llsraelian et al.l ( |2004|) sample. iBou vied (l2008l) has 
discussed the possibility of a different rotation history in EHS 
and NEHS to explain the lithium differences. In our models, ro- 
tation is treated in the same way for all the stars. In the third part, 
we discuss the influence of the presence of pla nets on the mixing 
in EHS , as suggested bv llsraelian et al.l(l2004l) and lChen & Zhaol 
(12006 ') and we present other models. Conclusions are given in 
the fourth section. 



2. Chemical peculiarities of exoplanet-host stars 

2.1. Overmetallicity of exoplanet-host stars 

The most important differe nce between NEHS and EHS is the 
overmetallicity of the later (lGonzaIezlll998t ISantos et af]l2000l 
2001, 2003). The mean overmetallicity of EHS is ~0.2 dex and 
the metallicity difference with NEHS is positive in 80% of the 
cases. 

Two scenarii have been proposed to explain this metal 
enrichment in central stars of planetary systems. The first 

scenario assumes a protostellar g as in itially metal -rich 

jPinsonneauIt et alj 120011; ISantos etalJl200lL l2003h . Following 
the "traditional" view, the gaseous giant planets are formed by 
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accretion of gas on planetesimals of around 10 earth masses. 
The higher the metalHcity (and the more dust particles there 
are), the more rapidly planetesimals are formed, and the higher 
is the probability to form giant planets before the dissipation 
of the protostellar disc. In this scenario, the star is overmetallic 
from the centre to the surface. 

The other suggestion is that the overmetallicity results from 
the accretion of planets, plan etesimals, asteroids , comets and 
dusts on the surface of the star (l Murravetal1l200lh . The mixing 
of matter occurring in the upper layers of the star homogenises 
the abundances, modifying the ratio between heavy elements, 
brought by the accreted objects, and hydrogen. If this mixing 
process was restricted to the convective zone, the metal enrich- 
ment would strongly depend on the stellar mass, w hich is not 
the case (Pinsonneault et al. 200ll [Santos et al.ll200ll) . However, 
IVauclaiii (2004) showed that, due to thermohaline mixing, the 
element dilution is more important in hotter stars than in cooler 
ones and could possibly account for a similar overabundance in 
stars of different masses. However, the possibility that the en- 
tire overmetallic matter sink inside the star is not negligible, in 
which case no overmetallicity would be left at all in the stellar 
outer layers. In any case, the large amount of accreted matter 
necessary to obtain the observed overmetallicity in cool stars, 
typically a hundred earth masses of iron, seems too large to be 
realistic. 

Another point is that the accretion scenario could only work 
if accretion still occurs while the star arrives at the beginning 
of the main sequence: the material accreted during the pre 
main-sequence phase is too much diluted in deep convective 
zones to lead to an observable effect. While the migration time 
scales are quite small, of order 10^ - 10^ years, the formation of 
giant plan ets can take a mu ch longer time, as large as a few 10^ 
years (e.g. llda & Linll2008h : planets may still form and migrate 
while the star is at the end of the PMS phase, but most of the 
accretion phase is expected to occur before that. 

Several studies have tried to differentiate the two sce- 

narii by differe nt manners, including asteroseismology 

jBazot & Vauclaiij |2004|) . A d eep analysis of th e star mu 
Arae, o bserved with HAR PS dBouchv et alj l2005ft has been 
given bv lBazot et alj (l2005l) . In spite of the excellent precision 
obtained for this sta r, it was not yet po ssible to disentangle 
the two possibilities. IVauclair et al.l (12008! ) gave evidence from 
asteroseismology that the exoplanet-host star t Hor, visible in 
the south hemisphere, has very probably been formed together 
with the Hyades cluster, inside the same overmetallic cloud. 
This gives further evidence that the accretion hypothesis as an 
explanation of the observed metallicity is ruled out. 

2.2. Lithium abundances in exoplanet-tiost stars 

In llsraehan et al.l (|2004|) . the authors present a comparison of the 
lithium abundances in EHS and NEHS. Two samples of stars are 
compared : a sa mple with 7 9 EHS and a comparison sample of 
157 NEHS fromChen etaD ( |2001|) . 

The lithium distributions in the samples show an impor- 
tant statistical excess of EHS having a lithium abundance 1.0 
< log e(Li) < 1.6, among the coolest stars. 

Lithium abundances of EHS with effective temperature be- 
tween 5850 and 6350 K are similar to those of the NEHS sample 
of|Chen et al., whereas at lower effective temperatures the EHS 
show a clear lithium underabundance relative to the comparison 
sample. The excess of lithium-poor EHS is concentrated in the 



range 5600 K < Tef f < 5850 K. This study was confirmed by 
IChen&Zhaol (12006 '). 

A different angular momentum history, due to the pres- 
ence of a more massive protostellar disc i n the case of 
EHS (lEdwards et alJ Il993t IWolffetalJ |2004 or to accre- 
tion of planetary-b odies into the star dSiess & Livid 119991 : 
llsraelian et al. !'2003), can affect the lithium abundances in EHS. 
A relative overmetallicity of the EHS, either initial or resulting 
from stellar pollution, could also be a clue to the abundance dif- 
ferences. 



3. Modelling 

The following study consists of comparing the evolution of 
lithium (a) in stella r models with a solar met allicity ("old" abun- 
dances as given bv lGrevesse & Noelslll993b . (b) in overmetallic 
models (with an initial overabundance of metals from the centre 
to the surface). The model with a solar metallicity is used as a 
reference model (standard model). 

3.1. Models and calibration 

The evolution of models with various m asses are computed 
with the T oulouse-Geneva Evolution Code dRichard et alj|1996t 
i Hui-Bon- Hoa 20Q2) and parameterized as follows : 

Inputs physics 

The OPAL2001 dRogers & Navfonovl |2002|) equation 
of state is implemented. We use the OPAL96 opaci- 
ties tables (Iglesias & RogersI 1 19961) completed by the 
[Alexander & FergussonI d 19941) low temperature opacities. 
For nuclear rea ctions rates, we use the analytical formulae 
of the NACRE (lAngulo et al.lll999l) compilation, taking into 
acc ount the three pp c h ains a nd the CNO tricycle, with the 
Bahcall & Pinsonneaulj dl992|) screening routine. Convection 
is treated according to the iBohm- Vitensel dl958h formalism 
of the mixing-length theory with ap = l/Hp - 1.75. For the 
atmosphere, we use a grey atmosphere following the Eddington 
relation. 

The abundance variations of the following chemical species 
are individually computed in the stellar evolution code: H, He, 
C, N, O, Ne and Mg. The heavier elements are gathered in Z . 
The initial composition follows the iGrevesse & Noelsl (Il993h 
mixture. 

Diffusion and rotation-induced mixing 

All models include gravita tional settling with d iffusion co- 
efficients computed as in iPaquette et al.l (Il986l) . Radiative 
accelerations are not computed here, as we only focus on 
solar-type stars where their effects are negligible. 



Rotation-induced mixing is computed as described in 
iTheado & Vauclairj d2003|). This prescription is an exten- 
sion of lZahnid 19921) and lMaeder & Zahnid 19981) . and introduces 
the feedback effect of the //-currents in the meridional circula- 
tion, due to the diffusion-induced molecular weig ht gradients. 
The e volution of the rotation profile follows the ISkumanichI 
(Il972h law. The models have a surface initial rotation velocity 
on the ZAMS equal to V,- - 100 km.s Other prescriptions, 
which include the angular mo mentum transport induced b y the 
mixing, have been g iven bv ICharbonnel & TalonI d 19991) and 
iPalacios et al.l (l2003l) . However, as rotation-induced mixing 
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alone cannot account for the flat rotation profile inside the Sun, 
these authors have later introduced the possible effect of internal 
grav ity waves triggered at the bottom of the convective zone 
(e.g. lTalon & Charbonnelll2005l) . Other authors suggest that the 
internal magnetic field is mor e important than internal w aves to 
transport angular momentum (iGough & Mclntvr3ll998l) . In any 
case, when adjusted on the solar case, all these prescriptions are 
able to reproduce the lithium depletion obs erved in the Hyades, 
and the results are finally qui te similar dTalon & Charbonnel 
fT998tlTheado & Vauclaidl2003[) . 

We also include a shear layer below the convective zone, treated 
as a tachocline (see Spiegel & Zahn 1992): this layer is param- 
eterized with an effective diffusion coeflic i ent de creasing expo- 
nentially downwards (see lBrun et al ] |1998lll999h : 

D,acho = £>f,cz exp ^In 2 ^ 

where Dbcz and rtcz are respectively the value of Dtadio at the 
bottom of the convective zone and the radius at this location, 
and A is a parameterized width. 
No overshooting is added here. 



3.2. Results 

3.2.1 . Influence of the metallicity on the lithium destruction 

Figures[T]presents the lithium destruction along the evolutionary 
tracks for overmetallic models, as a function of the effective tem- 
perature. A standard model (with a solar metallicity) is used as 
a reference. Here we compare the effect of metallicity (for both 
initial overmetallicity and accretion scenarii) on models with the 
same mass. Overmetallic models have a mass of 1.05 and 
initial metaflicities [Fe/H] = 0.12 0.18, 0.2 4, 0.30 and 0.36. The 
positions of the EHS studied in llsraelian et al. ( 2004,) are also 
plotted. For these stars (linear regression curve), the lithium de- 
pletion clearly increases for smaller effective temperatures. 

Stars with a higher surface metallicity have a deeper outer 
convective zone and lithium is more easily destroyed, due to the 
fact that the distance between the bottom of the convective zone 
and the nuclear destruction layer is smaller. 

The positions of the observed stars in Figure [1] shows that 
those which are cooler than the Sun can be accounted for by 
overmetallic models (this is also possible with other stellar 
masses and metallicities than showed here), whereas the stars 
hotter than the Sun cannot be accounted for in this framework. 

3.2.2. Influence of the stellar mass on the lithium destruction 

The lithium destruction in overmetallic models with the same 
surface metallicity but different masses is presented in Figure |2] 
as a function of effective temperature. The chosen metallicity is 
[Fe/H] = 0.24. As previously, a standard model is given as the 
reference. Overmetallic models have masses of 1.03, 1.04, 1.05, 
1.06 and 1.07 Mq 

Figure [3]presents the destruction of lithium in standard mod- 
els with masses of 0.98, 0.99, 1.00, 1.01 and 1.02 Mq. 

From Figure |2] we can check that the lower is the mass of 
the model, the more important is the lithium destruction during 
stellar evolution. 

As in Figure [T] the hottest observed stars cannot be ac- 
counted for by overmetallic models. 
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Fig. 1. Lithium destruction along the evolutionary tracks, as 
computed by TGEC, as a function of the effective temperature, 
for a standard model (solid line) of 1.00 M© and overmetallic 
models (dashed lines) of 1.05 M© and with initial metallicities 
[Fe/H] ^ 0.12 0.18 0.24, 0.30 and 0.36. The EHS observed by 
llsraehan et al.l (|2004|) are plotted (filled squares) and their metal- 
licities are indicated. Upper limits are filled triangles. The posi- 
tion of the Sun is indicated. The dot-dashed black line is the lin- 
ear regression of the lithium abundances for the observed stars. 
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Fig. 2. Same as Figure [T] but here overmetallic models with the 
same initial metallicity [Fe/H] - 0.24 and masses of 1.03, 1.04, 
1.05, 1.06 and 1.07 M© are compared. 



We can see from Figure [3] that the lithium destruction in 
these stars could be accounted for in standard models with 
masses between 1 .00 and 1 .02 M©, but this is inconsistent with 
their observed overmetallicities. 
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Fig. 3. Same as Figures [T] but for standard models (with solar 
metallicity) with masses of 0.98, 0.99, 1.00, 1.01 and 1.02 Mq. 



lithium destruction as observed in NEHS. For these reasons, we 
believe that the results that we obtain are significant. 



4. Extra mixing and planets migration 

Recent observations from I Chen & Zhaol (12006 ') confirm that 
the strong lithium destruction in EHS, observed for effective 
temperatures between 5600 and 5900 K does not result from 
the overmetallicity of these stars, but from the presence of 
the planets themselves. Indeed, the overmetallic NEHS do not 
present a strong lithium destruction, contrary to the EHS. As 
already evoked in Israelian et al. (2004): the strong lithium 
destruction in these stars could be associated with a late migra- 
tion of giant planets at the end of the pre-main sequence. As 
found from recent models ("e.g. Ilda & Linll2008l) . giant planets 
can still form at that time, and migrate towards the star by the 
so-called "type I" migration mechanism, which implies strong 
interaction with the disc. Angular momentum transfer from 
the planetary disc to the external outer layers of the star could 
induce differential rotation between the outer convective zone 
and the radiative zone below. This would lead to extra mixing in 
this region, during a short time scale, enough to destroy lithium 
more efficiently than in stars without planets. In the following, 
we show overmetallic models including such an extra mixing. 



3.3. Effective temperature scale 

The previous results depend strongly on the measures of the ef- 
fective temperatures of the EHS. If the effective temperatures 
were lower by 50 K, which is within the error bars, the hottest 
observed stars would be accounted for by the models. Several 
works compared the effective temperature determination using 
spectroscopy, as the sample presented here, w ith the determi- 
nation by the infrared flu x method (IRFM, s ee iBlackwell et al.l 
119801; lAlonso et al.lll996h . iRibas et all (l2003l) . using an IR pho- 
tometry method similar to IRFM, but without requiring the 
use of a bolometric calibration, found a excellent agreement 
in the entire temper a ture range w ith the spectroscopic anal- 
vsis of ISantosetal.1 (I2001L 
iRamf rez & Melendej d2004 



2003 ). The comparison done by 



2005h. using IRFM, seem s to im- 



ply th at the temperatures derived in lSantos et al.1 (l2004al) . as well 
as by iRibas et all (l2003h . are hotter by abou t 100 K than the 
tempe ratures they obtained. In a recent work, ICasagrande et al.l 
(|2006|) show that their IRFM ter nperatures agree very well with 
spectroscopic determi nations of lSantos et al.l (l2004ai l2005b and 
iLuck & Heite^ (l2005h . Thus, we can consider that the effective 
temperatures values used in this study have a good accuracy and 
that the possibility of systematic effects in spectroscopic metal- 
licity determinations is small. 

The uncertainties concerning the effective temperatures de- 
rived from our models are difficult to evaluate, in the same 
way as for all other 1-D, grey atmospheres computations. The 
TGEC has been tested together with six other stellar evolution 
codes in the Evolution and Seismic Tools Activity (ESTA) of 
the CoRoT mission (Monteiroet al. 2006). The differences for 
the effective temperature between the results given by all codes, 
using the same physics, is around one percent. In the present 
situation, using a more sophisticated atmosphere, or changing 
the internal physics could modify the derived effective tempera- 
ture. However these modifications would be similar for EHS and 
NEHS. The present models and their parameters have been ad- 
justed so as to precisely account for the seismic sun, and for the 



Figure m presents the profile of lithium destruction in over- 
metallic models of 1.05 and 1.07 Mq in different cases. The first 
case (continuous line) represents a star without planet, with a 
metallicity [Fe/H]o - 0.24 and an initial rotation velocity of 100 
km.s"', and the mixing described in Theado & Vauclair (2003*). 
This model is used as a "standard overmetallic model". 

The second case (dotted line) represents a similar model but 
with an initial rotation velocity of 200 km.s"'. We can see that 
the increase of the initial rotation velocity has a very weak in- 
fluence on the lithium depletion during the evolution. It is due 
to the efficient braking of the rotation, both models reaching the 
same rotation rate after 100 Myrs. 

In the third case (dashed line) we simulate the mixing asso- 
ciated with a shear instability below the convective zone, due to 
the angular momentum transfer, by an increase of the parame- 
ters of the shear layer in the standard overmetallic model during 
the beginning of the main sequence. According to ' Alibert et al.l 
(2005), migration timescales of giant planets around the central 
star are found to lie between 0.1 and 10 Myrs. Their simulation 
of formation process of one particular planet by accretion and 
migration around a Sun-like star show that the object evolves 
from a starting distance flj^,., = 15 AU to a final distance ~ 2.5 
AU from its star with a final mass of ~ 3.5 Mjup in 5.5 Myrs. In 
our model, the diffusion coefficient at the base of the convective 
zone is Dbc^ = 2.200x lO*" cm^.s"' (instead of D;,,.; = 2.600x10^ 
cm^.s ' in the standard overmetallic model) and the thickness of 
the shear layer is increased to A = 5.2 x 10^ cm (instead of 6.1 
X 10** cm for the standard overmetallic model), during the first 
4 Myrs, corresponding to the first 4 time steps of the evolution 
of the model. In Figure |4] the hottest observed stars can be ac- 
counted for by the model of 1.07 M© with an increase of the 
shear instability at the beginning of the main sequence. 

The parameters that we used to simulate this extra mixing 
layer have been chosen to obtain results consistent with the 
lithium depletion of the hotter observed stars. They are some- 
what arbitrary and certainly do not correspond to a unique so- 
lution. They only allow to show which kind of extra mixing is 
needed to explain the strong lithium depletion in cold EHS. A 
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complete hydrodynamical treatment of the planetary-migration- 
induced angular momentum transfer would be needed for a 
stronger conclusion. 
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Fig. 4. Same as Figures [T] to [3] but for standard overmetallic 
models ([Fe/H]o = 0.24 and V,- = 100 km.s"') in solid line, over- 
metallic models ([Fe/H]o = 0.24) with initial rotation velocity 
y, = 200 km.s"' in dotted line and overmetallic models ([Fe/H]o 
= 0.24) with an increase of the shear instability below the con- 
vective zone during the first 4 Myrs of the main sequence, in 
dashed line, with masses of 1.05 and 1.07 M©. 



instabilities inside the star. However, our results show that extra 
mixing in an initially overmetallic star is the best scenario to 
account for the extra-lithium depletion observed in EHS stars. 

Planetary migration is not the only possible explanation for 
the needed extra mixing in EHS. Another explanation of the ex- 
tra lithium depletion in EHS stars could be related to the rotation 
history of the s tar ( Edwards et a l. 1993; Wolff et al. 2004, 200^. 
More recently, iBouvierl (|2008|) has developed a consistent the- 
ory of the rotation history of solar type stars, taking into account 
three periods: the PMS, where the stellar rotation is magneti- 
cally coupled with the accretion disc, the approach to the ZAMS 
where the stellar rotation increases, and the Main Sequence re- 
laxation, where the stars spin down, leading to similar rotation 
rates for all the stars at the age of the Sun. The stellar rotation 
rate on the ZAMS depends on the initial rotation, but also on 
the lifetime of the disc: for larger lifetimes, the rotation rate is 
smaller To account for the observations of rotation in solar type 
stars from PMS to the sun, one has to assume that slowly rotating 
stars develop larger velocity gradients at the base of the convec- 
tive envelope than rapidly rotating stars. In this framework, the 
extra mixing needed to account for the larger lithium depletion 
in EHS than in NEHS could be related to a longer lifetime of the 
accretion disc, possibly needed for the formation and migration 
of giant planets. 

Such alternative explanations for the lithium observations in 
EHS versus non-EHS stars still present many unknowns and un- 
derlying assumptions. Although the planet migration scenario 
for the lithium extra depletion in EHS stars does not represent a 
unique possibility, the preliminary results we have obtained here 
are encouraging for further investigations. 
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5. Conclusion 

The purpose of this work was to present an analysis and interpre- 
tation of the observ ations of the lithium ab undances in EHS and 
NEHS as given by llsraelian et alj (|2004|) . We computed many 
evolutionary tracks with various masses and metallicities. 

The first result is that our models reproduce the increase of 
the lithium destruction below 5900 K towards smaller effective 
temperature. T he lithium depl e tion in stars cooler than the Sun, 
as observed by llsraelian et al.l ( |2004|) could be accounted for in 
the framework of overmetallic models. However, none of these 
models can account for the case of hotter stars: extra mixing 
below the convective zone is needed to account for the observa- 
tions. 

We have computed specific models including such an 
extra mixing. They show that the observed EHS hotter than 
the Sun can be accounted for in overmetallic models with an 
increase of the shearing instability below the convective zone 
during the first 4 Myrs of the main sequence. However these 
models are very simple and contain adjusted parameters. This 
resu lt is in favour of the suggestion by llsraelian et aLl (|2004|) 
and IChen & Zhaol (l2006h that the strong lithium depletion 
observed in these stars could be related to the planet migration 
mechanism, which would induce an extra mixing process below 
the outer convective zone at the end of the pre main-sequence. A 
complete treatment of this process would need hydrodynamical 
computations of the angular momentum transfer and shear flow 
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